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By monitoring multiple molecular transitions, force-clamp and trap-position-clamp methods have led to precise deter- 
minations of the free energies and free energy landscapes for molecular states populated in equilibrium at the same or 
similar forces. Here, we present a powerful new elaboration of the force-clamp and force-jump methods, applicable 
to transitions far from equilibrium. Specifically, we have implemented a live jump detection and force-clamp algo- 
rithm that intelligently adjusts and maintains the force on a single molecule in response to the measured state of that 
molecule. We are able to collect hundreds of individual molecular transitions at different forces, many times faster than 
previously, permitting us to accurately determine force-dependent lifetime distributions and reaction rates. Application 
of our method to unwinding and rewinding the nucleosome inner turn, using optical tweezers reveals experimental 
lifetime distributions that comprise a statistically-meaningful number of transitions, and that are accurately single ex- 
ponential. These measurements significantly reduce the error in the previously measured rates, and demonstrate the 
existence of a single, dominant free energy barrier at each force studied. A key benefit of the molecular yoyo method 
for nucleosomes is that it reduces as far as possible the time spent in the tangentially-bound state, which minimizes the 
loss of nucleosomes by dissociation. 



PACS numbers: 68.35.Ja, 61.25.Hq 

I. INTRODUCTION 

Single-molecule studies of biomolecules and their com- 
plexes have led to important new insights into the 
molecular mechanisms underlying many biological pro- 
cesses. Single-molecule force spectroscopy (SMFS), in 
particular, has proven a powerful method by which to 
study the mechanochemistry of molecular motors, includ- 
ing kinesinjnBI myosin J^^^ polymeraseJ^EI] helica seJ^^^, 
chromatin remodeller^SHIIl and the ribosomd^SMI. SMFS 
has also transformed our understanding of nucleic acid^^^ 
and protein folding^"* and biological assembly processes, 
such as chromatin compaction^^''^. In all of these cases, 
the key questions are: What microscopic states are populated 
at a given force? What are the characteristics of the force- 
dependent transitions among these states, including what is 
the distribution of lifetimes for each transition and what is 
the mean rate of each transition? These quantities determine 
the relevant free energy landscape of the biomolecules un- 
der study, and their diffusion constants for motion within this 
landscape. In practice, however, the level of detail resolvable 
in the free energy landscape may be limited by the number of 
molecular events that can feasibly be observed. 

To-date, the most detailed characterizations of biomolec- 
ular free energies and free energy landscapes have been 
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obtained in optical-tweezers-based force-clamfPHlSESEll 
and trap-position-clamfP32] experiments in which multiple 
molecular states are in equilibrium with each other at the 
same or nearly the same force. In such situations, it is 
often possible to observe hundreds of transitions back and 
forth among the different molecular states involved. Be- 
cause of the large number of transitions, the transition rates, 
and hence the free energies of the states, may be accurately 
determined^" '" ■ ^'^ By contrast, for molecular states 
that are not in equilibrium - that is, for molecular states 
for which spontaneous transitions back and forth are not ob- 
servable on experimental timescales at a single force - it is 
more challenging to collect a statistically large-enough data 
set to be able to determine the transition rates accurately. 
Most commonly under these circumstances, researchers carry 
out repeated measurements of the molecular force-versus- 
extension curve^'*'*"^*'^' in order to determine the distribu- 
tion of transition forces. To increase the rate of data acqui- 
sition in force-versus-extension measurements, tandem arrays 
of independent, identical molecules are often used^"* Nev- 
ertheless, each force-versus-extension curve contains a lim- 
ited number of transitions, and several practical limitations, 
such as surface sticking, tether rupture, protein dissociation, 
etc. mean that each molecule is only measurable for a finite 
period of time. Moreover, a fundamental limitation of the 
force-versus-extension approach, compared to force-clamp 
and trap-position-clamp methods, is that the force-dependent 
lifetime distributions and transition rates are not measured di- 
rectly. Instead, model-based approaches are required to de- 
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termine the transition rates from the measured distribution of 
transition forces, in which it is usually assumed that the tran- 
sition lifetimes are exponentially distributed at all force^SMI 
The force-jump method, in which a force-clamp alternates pe- 
riodically between a high-force and a low-force*"* over- 
comes a number of these difficulties. However, because the 
force-clamps must be maintained for long enough to ensure 
that the transitions have taken place, this method is not effi- 
cient. In addition, because the force must be maintained at 
a high value for an extended period of time, the force-jump 
method is especially vulnerable to tether rupture and other ir- 
reversible damage. 

The purpose of this paper is to present a powerful new elab- 
oration of the force-clamp and force-jump methods, which 
we call "the molecular yo-yo method", that permits the ef- 
ficient determination of lifetime distributions and transition 
rates for single molecules, even for molecular states that are 
not in equilibrium. The key innovation of the molecular yo-yo 
method to implement a live state-detection and force-clamp 
algorithm, that intelligently adjusts and maintains the force 
on a single molecule in real-time, in response to the measured 
state of the molecule. We describe two simple to implement 
live jump detection algorithms. The molecular yo-yo method 
is broadly applicable to out-of-equilibrium molecular transi- 
tions of all sorts. Here, to showcase the usefulness of this 
method, we report its implementation in experiments that seek 
an improved characterization of the unwinding and rewind- 
ing transitions of the nucleosome inner turn, which are out- 
of-equilibrium at near-physiological salt concentrationJSEl! 

The relevant states for nucleosome unwinding and rewind- 
ing are illustrated in Fig. In this paper, we probe the 
transition from state 1 to state as shown in Fig.[TJJ. Although 
states 1 and are out of equilibrium at any force, neverthe- 
less, using the molecular yo-yo method, we have been able to 
collect hundreds of molecular transitions, permitting us to ac- 
curately determine force-dependent lifetime distributions and 
force-dependent reaction rates for unwinding and rewinding 
the nucleosome inner turn. The molecular yoyo method is 
especially valuable for studies of nucleosomes because it re- 
duces as far as possible the time spent in the unwound state 
(state 0), correspondingly reducing the loss of nucleosomes 
by dissociation in a given period of time. Using this new 
method, we found that at each force studied, the lifetime dis- 
tributions are well-described as a single exponential, indicat- 
ing that the free energy landscapes relevant to winding and 
unwinding are each dominated by a single free energy bar- 
rier. We also demonstrate that the unwinding rates for tethers 
containing 2, 4 and 8 nucleosomes are accurately 2-, 4-, and 
8-fold faster, respectively, than for tethers containing a single 
nucleosome. This observation implies that nucleosomes on 
the same tether unwind independently, as has been previously 
assumed but not proven. Finally, we show that the rates mea- 
sured using the optical yo-yo method match those determined 
previously by the more laborious method of making a series of 
measurements at different forces, using a simple force clamp 
and force jumping. 
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FIG. 1. (A) Schematic of the four microstates for nucleosome 
unwinding/rewinding-'^. Histones are red. DNA is blue. For state 
2, the nucleosome is fully wrapped by nearly two turns of DNA. For 
state 1, the outer turn is unwrapped, but the inner turn is wrapped. 
For state 0, both the outer and inner turns are unwrapped, but the hi- 
stone octamer remains bound to the DNA. Finally, the histones may 
be unbound from the DNA. The shaded box encompasses the states 
and transitions probed in this paper. (B) Schematic illustration of the 
application of our molecular yo-yo method to the free energy land- 
scape for states not in equilibrium. The rate of transitions from state 
1 to state is non-negligible only at high force; the rate of transitions 
from state to state 1 is non-negligible only at low force. Follow- 
ing every transition, our algorithm detects the molecular state and 
adjusts the force so that the reverse transition can take place. By tog- 
gling between the high and low force configurations we are able to 
rapidly measure many transitions from a single molecule. 



II. MATERIALS AND METHODS 

A. DNA preparation and surface attachment 

Segments of DNA containing 4200 base-pairs, correspond- 
ing to a contour length of 1430 nm, were created from lin- 
earized pUC18 plasmids harboring an array of twelve 601 
nucleosome positioning sequences (NPSs), a kind gift from 
Dr. Daniela Rhoades^^ To facilitate robust attachments 
between the DNA and a microscope coverslip and between 
the DNA and the optically-trapped bead, the DNA is labeled 
with biotin on one end and an amine on the other end. To 
create DNA-tethered beads, the amine-labeled DNA end is 
covalently attached to the glass coverslip via a silane PEG 
N-hydroxysuccinamide linker while the biotin-labelled DNA 
end binds to a streptavidin coated polystyrene bead^^. 



B. Histone expression and purification 

pET vectors containing untagged Xenopus H2A, H2B, H3 
or H4 were a kind gift from Dr Karolin Luger. Histones 
are expressed in E. coli BL21-Gold(DE3), extracted from iso- 
lated inclusion bodies in buffer containing 7 M GuHCl, and 
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dialyzed into 8 M urea buffer. Histories are purified first by 
passage through a Q-sepharose column, and then bound to a 
Hi-Trap SP column (GE Healthcare Life Sciences), washed 
with 300 mM NaCl, and eluted with a step gi-adient to 600 mM 
NaCl. Finally, the histones are dialyzed into deionized dis- 
tilled H2O, lyophilized, and stored at — 80°C until needed. 

C. Nucleosome reconstitution 

Equimolar ratios of the four core histones are combined in 
buffer containing 7 M GuHCl and dialyzed into buffer con- 
taining 2 M NaCl, resulting in octamer formation^''. Af- 
ter isolation by gel filtration, octamer is mixed with carrier 
DNA - ultra pure salmon sperm DNA, sheared to 1000 bp 
(Invitrogen) - and continuously dialyzed into buffer without 
salt to form nucleosomes. We then assemble nucleosomes in 
situ by flowing a solution of nucleosomes, bound to carrier 
DNA, at 680 mM NaCl, into our optical tweezers flow cell, 
which consists of a flow channel cut out of double-sided sticky 
tape between a microscope coverslip and a microscope slide 
into which are drilled two holes for fluid inlet and outlet, re- 
spectively. In-situ nucleosomes exchange between the carrier 
DNA and the immobilized 601 DNA tether ensures occupancy 
of the 601 sites. We then flow in a 100 mM NaCl buffer, in 
preparation for optical tweezers measurements. For experi- 
ments discussed here, we prepared arrays of varying numbers 
of nucleosomes per tether in order to be able to compare rates 
of unwinding and rewinding with different number of nucleo- 
somes. 



D. Optical trapping instrumentation 

In the optical trapping setup used for these experiments, 
the beam from a 1064 nm laser (Ventus IR, Laser Quan- 
tum, Stockport, UK) is incident on an acousto-optic deflec- 
tor (AOD) (IntraAction DTD-274HA6), which serves to op- 
tically isolate the laser from the downstream optics, and to 
control the trapping laser power. Located between the AOD 
and the microscope objective (Nikon CFI xlOO, oil immer- 
sion, NA 1.25) is a telescope that expands the beam by a fac- 
tor of three to ensure that the back pupil of the objective is 
overfilled, as required for strong trapping. Beyond the ob- 
jective, the transmitted laser light is incident upon a quadrant 
photo-diode (QPD) (Phresh Photonics SiQu50-M), located in 
a plane conjugate to the back focal plane of the microscope 
condenser lens, where variations in the summed intensity of 
all four QPD quadrants are linearly proportional to the dis- 
placements of a trapped bead from the center of the trap along 
the beam direction. This method of determining^ bead posi- 
tion - "back-focal-plane interferometry" (BFPF/' - provides 
a sensitive measure of the force on the bead, which is propor- 
tional to the bead displacement from the trap center. 

The measurements described in this paper were carried 
out using an axial pulling geometry in which a piezo-electric 
translation stage moves the microscope coverslip along the di- 
rection of the laser beam, thus applying tension to a surface- 



tethered molecule in the axial direction and maintaining a sim- 
ple geometry at all extensions'^ 62 66 ^^^q ^gg Qf j-jjg 
ial pulling geometry, we have implemented a new calibration 
method that allows conversion from the experimental signals 
- stage displacement and scattering intensity - to calibrated 
values of the molecular extension and applied tension. Our 
calibration procedures are fully described in Ref. WT 

Our axial pulling geometry enables straightforward im- 
plementation of a reliable feedback loop that maintains a 
constant force on the tether^^. Specifically, the force mea- 
sured by the QPD is held constant by adjusting the posi- 
tion of a piezo-electric stage (NanoMAX311, Thorlabs). To 
achieve this force-clamp, we implemented a proportional- 
integral-derivative (PID) feedback controller using Lab View 
and Labview-MathScript, which carries out the conversion 
from QPD intensity and piezo-stage position to axial force and 
tether extension, followed by actuation of the piezo-stage PID 
control, at a cycle rate greater than 1000 Hz. Data acquisi- 
tion card PCIe-6343 (National Instruments) is used for both 
acquisition and output. 

When a nucleosome unwinds, the tether length increases by 
about 25 nm. Concomitantly, the force transiently decreases. 
However, within 10 ms, our force clamp has adjusted the trap 
position to increase the tension to the force-clamp value. This 
response time is limited by the mechanical response of the 
microscope stage, not by the computation time of our con- 
version algorithm. In fact, we can flexibly program an arbi- 
trary sequence of forces versus time, jumping out, if neces- 
sary, at programmable break points. This flexibility enables 
the molecular yo-yo method. 

In order to sustain meaningful measurements on the same 
molecular construct for extended periods of time, it is neces- 
sary to dynamically correct for any drift in the position of the 
piezo-stage and in the laser intensity. Since such drift is slow, 
it is satisfactory to apply a drift correction procedure every 
^200 s, which we do automatically. Specifically, stage drift 
is corrected for by measuring the position at which the bead 
contacts the coverslip, permitting us to update the piezo-stage 
calibration correspondingly. By holding the bead against the 
stage and measuring for 100 ms, we are able to establish the 
bead-coverslip separation to within about 2 nm. To correct for 
possible laser intensity drift, the force zero is established by 
placing the bead close to, but not in contact with, the cover- 
slip, so that the tether has a very low extension, and the cor- 
responding force is negligible. This procedure establishes the 
force to within 0.1 pN. 

III. RESULTS AND DISCUSSION 

A. Molecular yo-yo method 

The molecular yo-yo method operates as follows: At the 
start time, a high, unwinding force is suddenly applied to the 
nucleosome in the wound state (state 1). After a period of 
time at this force, the wound nucleosome transitions to the un- 
wound state (state 0), leading to an increase in extension. By 
using the preprogrammed force- versus -extension curve of the 
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FIG. 2. Unwinding and rewinding the nucleosome inner turn using tlie optical yo-yo method. (A) Extension versus time while a single 
nucleosome is repeatedly unwound at an unwinding force of 10.3 pN and rewound at a rewinding force of 3.8 pN over a period of about lis. 
Extension data obtained at the unwinding force is shown with a green background while the extension data obtained at the rewinding force is 
shown with an orange background. (B) Extension versus time, plotted over a restricted time range, showing four unwinding events from (A). 
(C) Extension versus time plotted for a further restricted time range, now showing a single molecular yo-yo cycle. The nucleosome unwinds 
at about 5.44 s and rewinds at about 6.02 s. 



states involved, the yo-yo algorithm recognizes the change in 
state, and, after a short delay, reduces the force to the rewind- 
ing force. After a period of time at the rewinding force, the 
nucleosome transitions back to the wound state (state 1), with 
a concomitant decrease in extension. The algorithm then rec- 
ognizes this state change and increases the force to the un- 
winding force once again. This cycle is then repeated multiple 
times, and at multiple unwinding and rewinding forces. Im- 
portantly in the nucleosome context, the molecular yo-yo al- 
gorithm reduces as far as possible histone-DNA dissociation 
from the unwound state by minimizing the time the nucleo- 
some spends in the unwound state. Each measured dwell time 
at a given force contributes to the lifetime distribution at that 
force. 

Fig. [2]\ shows a representative extension versus time trace 
that covers twenty-eight repeats of the molecular yo-yo live 
jump detection algorithm that detects and triggers on both 
unwinding and rewinding events. Fig. |2jj highlights nine of 
these cycles, and Fig.|2p a single cycle. In Fig.[2p, for times 
prior to 5.25 s, the nucleosome inner turn is wound and the 
force is 3.8 pN. At 5.25 s, the current yo-yo cycle is ini- 
tiated, when the force-clamp steps up to 10.3 pN, as indi- 
cated by the green background. This force jump is signaled 
by a corresponding stepwise increase in the extension trace 
of Fig. |2p, as the DNA tether stretches further in response 
to the increased force. The nucleosome remains wound at 



10.3 pN for approximately 0.2 s until it unwinds at 5.44 s, 
signaled by a jump in the extension of about 25 nm. Thus, 
during this particular cycle, the lifetime of the wound state at 
10.3 pN was measured to be 0.2 s. The nucleosome remains 
unwound at 10.3 pN for approximately 5 ms, which is the pe- 
riod of time required for the yo-yo algorithm to recognize the 
unwinding event. Once the transition to the unwound state has 
been recognized, the force-clamp steps down to 3.8 pN, as in- 
dicated by the orange background. The extension decreases 
concomitantly, because the DNA tether stretches less at the 
reduced force. The nucleosome remains unwound at 3.8 pN 
until 6.01 s before rewinding. This gives us a measurement of 
the lifetime of the unwound state at 3.8 pN to be 0.55 s. At 
6.02 s, rewinding is signaled by a decrease in the extension 
of about 22 nm. Once the rewinding transition is recognized, 
the next yo-yo cycle is then initiated by returning the force to 
10.3 pN. Triggering on both rewinding and unwinding events, 
as in Fig. [2] minimizes unnecessary measurements after the 
transition has occurred, and leads to a significantly faster data 
acquisition rate than otherwise would be possible. 

B. Live jump detection algorithm 

We have implemented two jump detection algorithms, one 
for unwinding and one for rewinding. Both jump detection 



5 




0.1 0.2 
probability 



I 1 .24 W^'^^y^^f'f'^'^^ff*'^ 



1.26 
1.24 
1.26 
1.24 
1.26 
1.24 
1.26 
1.24 
1.26 
1.24 
1.26 
1.24 
1.26 
1.24 
1.26 
1.24 




200 
time (ms) 



400 



FIG. 3. Extension versus time of individual nucleosomes lield at 
fixed force and their corresponding extension distributions. (A) Ex- 
tension versus time of a single nucleosome, field at a fixed unwinding 
force of 10.3 pN, measured every 1 ms and plotted as connected dots. 
At 3.6 s after the unwinding force is applied, an unwinding event is 
evident, in which the extension jumps by 25 nm. (B) Experimental 
probability distributions for the extension before (gray bars) and af- 
ter (white bars) the unwinding event. These histograms correspond 
to 3500 (gray) and 3300 (white) measurements of the extension with 
a bin size of 2 nm. The smooth curves are Gaussians, each plotted 
using the mean and standard deviation of the respective experimental 
distributions: Oi = 5.9 nm and ao = 5.6 nm, for state 1 and state 0, 
respectively, at this force. (C) Extension versus time of a single nu- 
cleosome, held at a fixed rewinding force of 3.8 pN, measured every 

1 ms and plotted as connected dots. At 0.98 s after the force is ap- 
plied, a rewinding event is observed as a stepwise 22 nm decrease in 
the extension. (D) Probability distributions determined before (gray 
bars) and after (white bars) the rewinding event with bins of size 

2 nm. These histograms correspond to 986 (gray) and 1314 (white) 
measurements of the extension. The smooth curves are Gaussians, 
each plotted using the mean and standard deviation of the respective 
experimental distributions: Cq = 5.9 nm and Oo = 6.2 nm, for state 
1 and state 0, respectively, at this force. 



algorithms rely on comparing live extension to preceding ex- 
tension measurements. Fig. [3]\ shows an example extension 
of a tether during a nucleosome unwinding event with no de- 
tection. At time zero, a force-clamp of 10.3 pN is applied 
to a single nucleosome. For a period of approximately 3.6 s, 
the nucleosome in state 1 remains stable at this force, while 
the extension fluctuates about a mean of about 1200 nm. At 
3.6 s, the nucleosome unwinds into state 0, signaled by an 
increase in extension caused by about 25 nm of DNA being 
released. Subsequently, state is stable, with its extension 
fluctuating about a mean of 1225 nm. hi Fig.|3jJ, histograms 



FIG. 4. Examples of individual nucleosome unwinding event detec- 
tion. (A) Eight extension versus time traces at an unwinding force of 
10.3 pN, each showing an example unwinding event recognized by 
our jump detection method. In each panel, the black dots connected 
by gray lines represent the extension measurement. The lower dotted 
line corresponds to the mean of the ni (40) points measured before 
the jump, and represents the extension of state 1, the wound state. 
The upper dotted line represents the extension of state 0, the un- 
wound state, taken as the wound states extension plus 25 nm. The 
dashed line is 10 nm from the lower dotted line (state 1), and repre- 
sents the threshold extension for jump detection in each trace. (B) 
Count of the number of unwinding events versus time, determined 
from the eight example traces shown in (A). 



of the measured extension, represented as probabilities, before 
(gray) and after (white) the unwinding transition agree well 
with overlaid Gaussian distributions, plotted using the mean 
and standard deviations of the respective measured extension 
distributions. 

Similarly, Fig. [3p shows an example nucleosome rewind- 
ing event with no detection. At time zero, a force-clamp of 
3.8 pN is applied to a single nucleosome, initially in state 
(unwound). For a period of approximately 0.98 s, the nucleo- 
some in state remains stable at this force, while the extension 
fluctuates about a mean of about 1163 nm. At approximately 
0.98 s, the nucleosome rewinds into state 1, signaled by a de- 
crease in extension of about 22 nm. Subsequently, state 1 is 
stable, with its extension fluctuating about a mean of 1 141 nm. 
In Fig.[3p, histograms of the measured extension, represented 
as probabilities, before (gray) and after (white) the unwind- 
ing transition agree well with overlaid Gaussian distributions, 
plotted using the mean and standard deviations of the respec- 
tive measured extension distributions. These data were col- 
lected at a rate of 1 kHz, where extension fluctuations effec- 



lively correspond to uncorrelated, white noise. 

To determine in real-time whether an unwinding transition 
has occurred, our yo-yo algorithm looks back at the previous 
«i +«2 extension measurements. If all of the immediately 
preceding ni measurements exceed the mean of the 112 mea- 
surements previous to those ni measurements by more than 
a threshold value (Ai), then the yo-yo algorithm recognizes 
that a jump in extension has occurred. To assess this simple 
scheme for false positives, we inquire: What is the probability 
that the algorithm recognizes a transition from state 1 to state 
0, when in fact none occurred? The probability that a single 
point exceeds the threshold, when the nucleosome remains in 
state 1, may readily be seen to be equal to 
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where erfc is the complementary error function, Oy is the stan- 
dard deviation of the extension fluctuations in state 1, and Ai 
is the difference in extension between the threshold extension 
and the mean extension of state 1. The factor (1 + l/«2) ac- 
counts for the expected variance in the mean extension de- 
termined from n2 measurements. All of the unwinding yo- 
yo measurements presented in this paper employed n\ = 5 
and ¥12 — 40. With <3\ =5.9 nm, as found experimentally 
at 10.3 pN, it follows that P = 0.031. Given flie specified 
conditions for recognizing the transition, the probability of a 
false positive is P"', leading to a false negative probability of 
(0.031)5 =2.9xl0-^ which we consider entirely acceptable. 
Similar considerations apply to unwinding at other forces. 

Fig. |4|\ shows eight example nucleosome unwinding 
events, as detected by our jump detection algorithm. In each 
case, the measurement begins with the nucleosome in the state 
1 (wound) and with the application of the unwinding force of 
10.3 pN at time zero. After a variable time interval, the nu- 
cleosome unwinds into state 0. In each case, after about 5 ms, 
corresponding to ni = 5, the jump detection algorithm reg- 
isters the unwinding transition and reduces the force, having 
measured the lifetime of state 1 at the unwinding force. In all 
cases, the transition is abrupt and is promptly detected by the 
algorithm. In each panel, the lower dotted line corresponds 
to the mean of the previous n2 = 40 points before the jump, 
and represents the extension corresponding to state 1 . The up- 
per dotted line corresponds to the mean extension of the un- 
wound state at an extension 25 nm higher than the lower dot- 
ted line. The dashed line represents the threshold, described 
previously, and is located Ai = 11 nm above the mean exten- 
sion of state 1. Fig.|4^ summarizes the cumulative number of 
unwinding events versus time, determined from the measure- 
ments of Fig. |4jA. 

We use a modified jump detection algorithm for nucleo- 
some rewinding. At the beginning of this protocol, the exten- 
sion of a nucleosome in state 1 is determined by a 1 s exten- 
sion measurement at the rewinding force. The extension of 
state at the rewinding force is then taken to be the extension 
of state 1 plus 23 nm. Once the force is jumped to the rewind- 
ing force, our algorithm registers a rewinding event when the 
mean of the previous rii extension measurements lies within a 
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FIG. 5. Examples of individual nucleosome rewinding event de- 
tection. (A) Eight extension versus time traces at a rewinding force 
of 3.8 pN, each showing an example rewinding event recognized by 
our rewinding jump detection method. In each panel, the black dots 
connected by gray lines represent the extension measurement. The 
lower dotted line corresponds to the previously-determined extension 
of state 1 at 3.8 pN. The upper dotted line represents the extension of 
state 0, taken as the extension of state 1 plus 23 nm. The dashed line, 
10 nm from the lower dotted line (state 1), represents the threshold 
extension for jump detection in these data. (B) Count of the number 
of rewinding events versus time, determined from the eight example 
traces shown in (A). 



threshold value, Ai, of the extension of state 1. In this case, 
the probability of falsely identifying a transition, when none 
has occuiTed, is 



1 ( ti.\Jn{ 
P = -erfc 
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Fig. |5}\ shows eight example nucleosome rewinding events, 
as detected by this algorithm. In each case, the measurement 
begins with the nucleosome in the state and with the ap- 
plication of the rewinding force of 3.8 pN at time zero. Af- 
ter a variable time interval, the nucleosome rewinds into state 
1. Then, after a time that can be seen to vary from 10 to 40 
ms, the jump detection algorithm registers the rewinding tran- 
sition and increases the force, having measured the lifetime 
of state at the rewinding force. In all cases, the transition 
is unambiguous and is readily detected by the algorithm, al- 
though not as rapidly as in the case of unwinding. For these 
example rewinding events, for which a = 6 nm, n\ = 40, and 
Ai = 10 nm, a false-positive jump detection occurs with a 
probability of 3 x 10^^^, effectively impossible on experimen- 
tal time-scales. Fig.lSfe summarizes the cumulative number of 
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time (s) 

FIG. 6. Cumulative fraction of nucleosomes unwound versus time at 
10.3 pN determined on tethers with 8 (green), 4 (blue), 2 (red), or 
1 (purple) nucleosomes. Each curve corresponds to measurements 
performed on a single tether. Each step in the staircases corresponds 
to an independent unwinding event. The smooth curves are exponen- 
tials with rates Sk, 4k, 2k, and k, where k is determined by the max- 
imum likelihood method from the single nucleosome dat^^l. The 
apparent rate scales with the number of nucleosomes. Measurements 
on different tethers give the same rate for single nucleosome unwind- 
ing. 

rewinding versus time, determined from the measurements of 
Fig.|5]\. 

C. Distribution of lifetimes at fixed force 

To determine the Hfetime distribution of state 1 or state 
at each force studied, we count the cumulative number of un- 
winding or rewinding events as a function of time, as indicated 
in Fig. |4j3 and Fig. |5j3. Normalizing by the total number of 
counts, these collections of unwinding and rewinding times 
determine the cumulative fraction of nucleosome inner turns 
unwound and wound, respectively, as a function of time, at a 
given force. In Fig.[6j we show unwinding data at 10.3 pN ob- 
tained on tethers with one, two, four and eight nucleosomes in 
purple, red, blue and green respectively. We plot two curves 
each for 1, 2, 4, and 8 nucleosomes. Each was obtained in 
separate experiments with different tethers and nucleosomes, 
demonstrating individual tether-to-tether and nucleosome-to- 
nucleosome repeatability. 

If we assume first-order unwinding (rewinding) kinetics, 
the probability that a nucleosome initially in state 1 (or 0) will 
have undergone unwinding to state (or 1) within a time t 
after initiation of the force clamp is given by an exponential 
function; 

p^l~e-'\ (3) 

where k is the unwinding (rewinding) rate. For independent 
nucleosomes, we expect the apparent unwinding rate for 2, 4, 
and 8 nucleosomes to be 2 times, 4 times and 8 times, respec- 
tively, the rate for a single nucleosome. The smooth curves 




time (s) 



FIG. 7. Cumulative fraction of nucleosomes unwound versus time 
for unwinding and rewinding determined using the molecular yo-yo 
method. (A) Distribution of unwinding times at 8.5 pN (red), 9.4 pN 
(green), and 10.3 pN (blue) shown using a logarithmic time axis, with 
disitrubutions determined using 201, 211, and 354 unwinding events, 
respectively. Each step is an independent unwinding event. (B) Dis- 
tributions of rewinding times at 3.3 pN (orange) and 3.8 pN (purple), 
after unwinding at 14. 1 pN, determined using 66 and 78 rewinding 
events, respectively. In both (A) and (B), the corresponding, over- 
layed smooth black curves are exponentials with rates determined by 
the maximum likelihood method. 

in Fig. |6] correspond to exponentials with rate set according 
to this rule, using the maximum likelihood value of the un- 
winding rate determined from data obtained on a single nu- 
cleosome. Evidently, the model curves provide an excellent 
description of the behavior with 1, 2, 4, and 8 nucleosomes. 
This observation directly demonstrates that at the forces stud- 
ied nucleosomes on the same DNA tether unwind indepen- 
dently, as has previously been assumed but not proven. 

Fig. [TjA. displays the distributions of unwinding times, us- 
ing a logarithmic time axis, at 8.5, 9.4, and 10.3 pN, each 
determined from multiple nucleosomes, collected together to 
yield a single distribution at each force. Similarly, Fig. |7jj 
displays the distribution of rewinding times at 3.3 and 3.8 pN 
for nucleosomes, that were unwound at 14.1 pN. For both 
Fig. [7|\ and B, each distribution is represented as the frac- 
tion of nucleosomes unwound. The collected events at each 
force constitute a sufficiently large data set (201, 211, 354, 
66 and 78 transitions at 8.5, 9.4, 10.3, 3.3 and 3.8 pN, re- 
spectively) to enable us to not only determine the unwinding 
rate at the force in question but also to test whether an expo- 
nential distribution of lifetimes is a correct description. The 
solid black lines in Fig.|7}\ and B correspond to EQ.|3]calcu- 
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lated using the maximum likelihood values of the unwinding 
and rewinding rates respectively. Evidently, this model pro- 
vides an excellent description of our measured lifetime distri- 
butions with zero adjustable parameters, indicating that a sin- 
gle exponential lifetime distribution is the correct description, 
and that the transition rates are 0.062 ±0.004, 0.16 ±0.01, 
0.40±0.02, 1.2±0.1, and 0.38±0.04 i"' at 8.5, 9.4, 10.3, 
3.3, and 3.8 pN, respectively. To objectively assess how well 
EQ. |3] accounts for the measured lifetime distributions, we 
have binned the unwinding lifetime measurements shown in 
Fig. [t] into logarithmically-sized bins^^. Binning the data to 
obtain the distribution of lifetimes ensures that the number of 
counts in the different bins are statistically independent of one 
another, which is not the case for the cumulative distributions 
of Figs. |6] and |7] It also permits us to simply determine the 
standard error for each bin as the square-root of the number of 
counts in each bin. The corresponding histograms, including 
error bars, are compared with the model distribution corre- 
sponding to EQ.[3]in Fig. [8] To determine the goodness of fit, 
we calculated the reduced chi-squared: 
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where the sum runs from 1 to n = 13 bins, O, is the observed 
number of counts in bin /, and Ei is the expected number of 
counts in bin W^. The reduced x^-values are 0.58, 1.3, and 
1.21, for the unwinding distributions at 8.5, 9.4, and 10.3 pN, 
respectively. For the rewinding distribution at 3.3 and 3.8 pN, 
the reduced x^-values are 1.0 and 1.3, respectively. Thus, in 
every case, is close to unity, indicating that there is not 
a statistically significant deviation between the data and the 
model of EQ.[3]at any of the forces studied. 



D. Force dependent unwinding and rewinding rates of 
the nucleosome inner turn 



FIG. 8. Histogram of nucleosome unwinding times shown using a 
logarithmic time axis. Unwinding times at 8.5 pN (red), 9.4 pN 
(blue), and 10.3 pN (green) were logarithmically-binned with dis- 
tributions deteiTnined using 201, 21 1, and 354 unwinding events, re- 
spectively. Each bin is a factor of 3 longer than the preceding bin. 
The histogram of the number of events with lifetimes in each time 
bin is displayed as a staircase. The standard error for the number of 
counts in each bin, shown as the error bar, is determined by count- 
ing (Poisson) statistics. An exponential distribution, corresponding 
to the characteristic unwinding rate, determined by the maximum 
likelihood method, is overlaid as the thin continuous line. 



Fig. |9] summarizes our measurements of the force- 
dependent rates of unwinding and rewinding the nucleosome 
inner turn, obtained with the molecular yo-yo method. For 
comparison. Fig. |9] also includes our previously published 
results for these rateJS! Evidently, measurements of the 
nucleosome inner turn unwinding and rewinding rates, ob- 
tained using the molecular yo-yo method, show good agree- 
ment with those obtained previously. At a given force, how- 
ever, the number of transitions available via the yo-yo method 
is several-fold larger than the number previously available, 
yielding values for the transition rates that are accurate to 
within about 10%. Thus, in future nucleosome studies the 
molecular yo-yo method will enable us to resolve subtle dif- 
ferences in the kinetics of different nucleosomes, with differ- 
ent histone variants and modifications, for example. 

E. Conclusion 

In this paper, we have introduced and demonstrated a pow- 
erful elaboration of the force clamp method, called the molec- 



ular yo-yo, that is broadly applicable to molecular transitions 
that are far from equilibrium. The molecular yo-yo method 
implements a live jump-detection and force-clamp algorithm, 
that intelligently adjusts and maintains the force on a single 
molecule, in response to the measured state of that molecule. 
Thus, we are able to realize hundreds of individual molecular 
transitions between molecular states at different forces, per- 
mitting us to accurately determine force-dependent lifetime 
distributions and reaction rates. Compared to force-versus- 
extension measurements, the molecular yo-yo method directly 
measures these key quantities, while maximizing data acqui- 
sition rate and efficiency. Compared to force-jump measure- 
ments, the molecular yo-yo method minimizes the time spent 
at high force, where molecular complexes can dissociate. We 
presented measurements detailing the specific application of 
the molecular yo-yo to unwinding and rewinding the nucleo- 
some inner turn, using optical tweezers. Because the molec- 
ular yo-yo minimizes the time spent in the unwound state 
(state 0), nucleosome dissociation is also minimized, permit- 
ting hundreds of transitions to be obtained from a single con- 
struct before nucleosome dissociation ends the measurement. 
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FIG. 9. Unwinding and rewinding rates of the nucleosome inner turn 
plotted as a function of force. Sliown as blue circles are the unwind- 
ing rates obtained from 12 nucleosomes undergoing a total of 1086 
unwinding events with the molecular yo-yo method. Shown as blue 
triangles are the rewinding rates, obtained using the molecular yo-yo 
method, from 6 nucleosomes undergoing a total of 167 rewinding 
events. The previously published unwinding and rewinding rates are 
shown as red circles and triangles, respectively. These data were ac- 
quired using a total of 334 nucleosomes with 265 unwinding events 
and 136 rewinding events. 



Our molecular yo-yo measurements of unwinding and rewind- 
ing the nucleosome inner turn reveal experimental lifetime 
distributions that are accurately single exponential, indicating 
the existence of a single dominant free energy barrier between 
states 1 and 0. We also demonstrate that the unwinding rates 
for tethers containing 2, 4, and 8 nucleosomes are accurately 
2-, 4-, and 8-fold faster, respectively, than for tethers contain- 
ing a single nucleosome. This observation implies that nucle- 
osomes on the same tether unwind independently, as has been 
previously assumed but not proven. Finally, we note that im- 
proved throughput is not the only benefit allotted by live jump 
detection, only the most obvious. A variety of more elaborate 
experimental force protocols enabling new measurements are 
clearly made possible with variations of the technique pre- 
sented here. Additionally, live jump detection methods are 
transferable to other SMFS methods, such as those that em- 
ploy an atomic force microscope. 
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